Observation of the scissors mode and superfluidity of a trapped Bose-Einstein 

condensed gas 
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We report the observation of the scissors mode of a Bose- 
Einstein condensed gas of 87 Rb atoms in a magnetic trap, 
which gives direct evidence of superfluidity in this system. 
The scissors mode of oscillation is excited by a sudden rotation 
of the anisotropic trapping potential. For a gas above T c 
(normal fluid) we detect the occurrence of oscillations at two 
frequencies, with the lower frequency corresponding to the 
rigid body value of the moment of inertia. Well below T c the 
condensate oscillates at a single frequency, without damping, 
as expected for a superfluid. 

PACS numbers: 03.75.Fi, 05.30. Jp, 32.80.Pj, 67.90.4-z 



The relationship between Bose-Einstein condensation 
(BEC) [0 and superfluidity has been studied extensively 
in liquid helium 0, but only recently has it been pos- 
sible to examine it in condensates of dilute alkali metal 
vapours ||,|| . Helium below its critical point is described 
by a two- fluid model and the liquid is endowed with a new 
degree of freedom, namely the relative motion between 
the normal fluid and the superfluid. "This degree of free- 
dom is the essence of the transport phenomena in He II 
known collectively as superfluidity" H . In the dilute al- 
kali metal condensates various phenomena which imply 
the occurrence of superfluidity have been observed e.g. 
collective modes of excitation || and demonstrations of 
the coherence of the wavefunction |Q. However in a re- 
cent theoretical paper D. Guery-Odelin and S. Stringari 
H describe how the superfluidity of a trapped BEC may 
be demonstrated directly and we report the results of 
such an experiment. Guery-Odelin and Stringari analyse 
the so-called scissors mode in which the atomic cloud os- 
cillates with respect to the symmetry axis of the confining 
potential and they point out that the scissors mode has 
been used in nuclear physics to demonstrate the super- 
fluidity of neutron and proton clouds in deformed nuclei 

IS 

The full theoretical analysis of the scissors mode is 
given in || and we only outline the key points here. The 
starting point is a BEC in an anisotropic harmonic po- 
tential with three different frequencies lu x « u> y < oj z . 
The scissors mode may be initiated by a sudden rotation 
of the trapping potential through a small angle as indi- 
cated in Fig.|l]. In the subsequent motion, the cloud is 
not deformed provided that the change in the potential 



is too small to excite shape oscillations. For a thermal 
gas both rotational and irrotational fluid flow occur in 
the scissors mode and the normal fluid is predicted to 
exhibit two frequencies corresponding to those forms of 
motion. For the BEC there is only irrotational flow be- 
cause of its single valued wavefunction and therefore it 
only exhibits one frequency, which is different from cither 
of the frequencies observed for the thermal cloud. 

In our experiment the trapping potential is created by 
a time-averaged orbiting potential (TOP) trap jll] which 
is a combination of a static quadrupole field, of gradient 
Bq in the radial direction, and a time-varying field 

Br (t) = B r (cos fit e x + sin fit e y ) + B z cos fit e z . 

The term B z (t) — B z cos fit e z is additional to the usual 
field of amplitude B r rotating in the xy-plane. The effect 
of the additional term B z (t) is to tilt the plane of the 
locus of B = by an angle £ = tan~ 1 (B z / B r ) with 
respect to the xy-plane. This causes the symmetry axes 
of the potential to rotate through an angle <p w |£ in 
the xz-plane (this analytic result is only valid for £ 2 <C 
1). It is only necessary to consider the xz-plane since 
the absorption image of the cloud is projected onto this 
plane. Tilting the locus of B = reduces the oscillation 
frequency in the z direction from its value when B z = 0. 
Thus simply switching on B 2 (t) also changes the cloud 
shape and so excites quadrupole mode oscillations. To 
avoid this we first adiabatically modify the usual TOP 
trap to a tilted trap and then quickly change B z (£) to 
— B 2 (t) . This procedure rotates the symmetry axes of the 
trap potential by 2<f> without affecting the trap oscillation 
frequencies (Fig.|l|). 

Our apparatus for producing BEC of 87 Rb is described 
in p2[ and only a few relevant features are outlined be- 
low. We trap and cool the atoms in a small glass cell 
that forms part of a differentially pumped system with 
two magneto-optical traps (MOT) E|, one of which is 
a pyramidal configuration of mirrors inside the vacuum 
system [M. The Helmholtz coils which create the oscil- 
lating fields along the x, y and z directions are driven by 
audio amplifiers and it was found to be extremely im- 
portant to filter out high frequency noise, especially in 
the coils driving the field along the z-axis jl5j . The 87 Rb 
atoms in the F = 2, mp = 2 state, were probed by a 10 /is 
pulse of laser light propagating along the y-axis, resonant 
with the F = 2,mp = 2 — > F' = 3, m' F = 3 transition. 
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The probe pulse was synchronized to the point in the 
rotation when the magnetic field points along the y-axis 
and this probing scheme is not affected by B z (t). We 
are restricted to probing at times which are multiples of 
the rotation period of the field (143 /is) but this is very 
much smaller than the period of the oscillations in the 
trap, as it must be for the TOP trap to work. 

The following experimental procedure was used to ex- 
cite the scissors mode both in the thermal cloud and in 
the BEC. Laser cooled atoms were loaded into the mag- 
netic trap and after evaporative cooling the trap frequen- 
cies were u> x = 90 ± 0.2 Hz and lu z = y/8 uj x . The trap 
was then adiabatically tilted in 1 s by linearly ramping 
B z (t) to its final value, corresponding to 4> = 3.6° and a 
reduction of the trap frequency uj z by ~ 2%. Suddenly 
reversing the sign of B 2 (t) in less than 100 /is excites 
the scissors mode in the trapping potential with its sym- 
metry axes now tilted by — </>, as described above. The 
initial orientation of the cloud with respect to the new 
axis is #o = 2(/>, so this angle is the expected amplitude 
of the oscillations (Fig||). The angle of the cloud was 
extracted from a 2-dimensional Gaussian fit of the ab- 
sorption profiles such as those shown in Fig.0. 

For the observation of the thermal cloud the atoms 
were evaporatively cooled to 1 /xK which is about 5 times 
T c , the temperature at which quantum degeneracy is ob- 
served. At this stage there were ~ 10 5 atoms remaining 
with a peak density of no ~ 2 • 10 12 cm~ 3 . The scissors 
mode was then excited and pictures of the atom cloud in 
the trap were taken after a variable delay. The results of 
many runs are presented in Fig.^|(a) showing the way the 
thermal cloud angle changes with time. The model used 
to fit this evolution is the sum of two cosines, oscillat- 
ing at frequencies u>i and ll>2- From the data we deduce 
wi/2tt = 338.5 ± 0.8 Hz and w 2 /2tt = 159.1 ± 0.8 Hz. 
These values are in very good agreement with the val- 
ues 339 ± 3 Hz and 159 ± 2 Hz predicted by theory ||; 
which correspond to lo\ = lu z +lu x and L02 — uj z — u x . We 
measured uj x and lo z by observing the center of mass os- 
cillations of a thermal cloud in the untilted TOP trap and 
calculating the modification of these frequencies caused 
by the tilt. The amplitudes of the two cosines were found 
to be the same, showing that the energy is shared equally 
between the two modes of oscillation. 

To observe the scissors mode in a Bose-Einstein con- 
densed gas, we carry out the full evaporative cooling 
ramp to well below the critical temperature, where no 
thermal cloud component is observable, leaving more 
than 10 atoms in a pure condensate. After exciting the 
scissors mode we allow the BEC to evolve in the trap for 
a variable time and then use the time-of-flight technique 
to image the condensate 15 ms after releasing it from 
the trap. The repulsive mean-field interactions cause the 
cloud to expand rapidly when the confining potential is 
switched off, so that its spread is much greater than the 
initial size. The aspect ratio of the expanded cloud is op- 



posite to that of the original condensate in the trap, so 
that the long axis of the time-of-flight distribution is at 
90° to that of the thermal cloud as shown in FigJ|. How- 
ever this difference in the orientation does not affect the 
amplitude of the angle of oscillation. The scissors mode 
in the condensate is described by an angle oscillating at 
a single frequency uj c 

9 (t) = -<f> + 9 cos (uj c t) (1) 

Figure || (b) shows some of the data obtained by exciting 
the scissors mode in the condensate. Consistent data, 
showing no damping, was recorded for times up to 100 
ms. From an optimized fit to all the data for the function 
in Eq.[l] we find a frequency of lu c /2tt = 265.6 ± 0.8 Hz 
which agrees very well with the predicted frequency of 
265 ± 2 Hz from uj c = y/u) x + w\. The aspect ratio of 
the time-of-flight distribution is constant throughout the 
data run confirming that there are no shape oscillations 
and that the initial velocity of a condensate (proportional 
to 9) does not have a significant effect. 

These observations of the scissors mode clearly demon- 
strate the superfluidity of Bose-Einstein condensed ru- 
bidium atoms in the way predicted by Guery-Odelin and 
Stringari || . Direct comparison of the thermal cloud and 
BEC under the same trapping conditions shows a clear 
difference in behaviour between the irrotational quan- 
tum fluid and a classical gas. Another distinction is the 
lack of damping in the superfluid. The damping of the 
classical gas is not apparent from our data because it 
is not sufficiently dense and the time between collisions 
is many oscillation periods. However the damping in a 
thermal cloud can be calculated using the Direct Simu- 
lation Monte Carlo method |l6j, for a mean density and 
temperature that are roughly the same as those of the 
BEC. (Since the condensate and the thermal cloud have 
different spatial distributions the comparison is only ap- 
proximate.) The results of such a numerical calculation 
for a thermal cloud at a temperature of 90 nK and a 
peak density of 2 • 10 14 cm~ 3 are shown in FigJJ and 
the plot shows that the scissors mode of a classical gas 
is strongly damped under these conditions. Note that 
even for such high densities both frequency components 
still occur in the thermal cloud, giving behaviour that 
is clearly different from the undamped, single frequency 
oscillation observed for the BEC. The results of this nu- 
merical simulation show that the amplitude of the lower 
frequency component is much smaller than that of the 
higher frequency one. The higher frequency tends to- 
wards oj c , the same frequency as the condensate, when 
the density increases so that the hydrodynamic regime is 
reached (where there are many collisions per oscillation 
period). However in this regime the damping is so strong 
that only a few oscillation periods would be observed as 
shown in Q . 

In the near future we plan to measure the frequency 
of the scissors mode at finite temperatures in the con- 
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densate i.e. where a trapped thermal cloud is present 
in addition to the BEC. Under these conditions the scis- 
sors mode should be damped S in a similar way to the 
quadrupole oscillations at finite temperature [|Tt|] . There 
is still not good agreement between the observed change 
in frequency of the m=0 oscillation mode of a conden- 
sate at finite temperature and theoretical predictions 
and therefore using the scissors mode as an alternative 
means of probing the interaction between the condensate 
and thermal cloud is important. We have found that the 
measurement of the angle of orientation of the cloud in 
the scissors mode is quite robust and could be applied in 
situations close to T c where only a small proportion of 
the atoms are condensed. In this way the scissors mode 
should enable studies of the relative motion of the ther- 
mal cloud and condensate corresponding to motion of the 
superfluid through a normal fluid. 
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program (No. ERB FMRX-CT96-0002). O.M. Marago 
acknowledges the support of a Marie Curie Fellowship, 
TMR program (No. ERB FMBI-CT98-3077). 




FIG. 1. The method of exciting the scissors mode by a sud- 
den rotation of the trapping potential. The solid lines indicate 
the shape of the atomic cloud and its major axes. The dotted 
lines indicate the shape of the potential and its major axes, 
(a) The initial situation after adiabatically ramping on the 
field in z direction, with cloud and potential aligned, (b) The 
configuration immediately after rotating the potential, with 
the cloud displaced from its equilibrium position, (c) The 
large arrow indicates the direction of the scissors mode oscil- 
lation and the smaller arrows show the expected quadrupolar 
flow pattern in the case of a BEC. The cloud is in the middle 
of an oscillation period. (The angles have been exaggerated for 
clarity.) 





100 um 

FIG. 2. (a) Image of the thermal cloud in the trap, (b) 
Image of the condensate after 15 ms of time of flight. The 
condensate expands most rapidly along the direction which is 
initially most tightly confined, leading to a 90° difference in 
orientation as compared to the thermal cloud. 
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FIG. 3. (a) The evolution of the scissors mode oscillation 
with time for a thermal cloud. For a classical gas the scissors 
mode is characterized by two frequencies of oscillation. The 
temperature and density of our thermal cloud are such that 
there are few collisions, so no damping of the oscillations is 
visible, (b) The evolution of the scissors mode oscillation for 
the condensate on the same time scales as the data in (a). 
For the BEC there is an undamped oscillation at a single 
frequency ui c . This frequency is not the same as either of the 
thermal cloud frequencies. 
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FIG. 4. A numerical simulation of the scissors mode os- 
cillation for a thermal cloud with a temperature and density 
comparable to the condensate in our experiment. The two 
frequency components are present and the damping time is 
about 15 ms. 
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